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ABSTRACT 

Eollowing our previous work, which related generic features in the sky-averaged (global) 21-cm signal to 
properties of the intergalactic medium, we now investigate the prospects for constraining a simple galaxy 
formation model with current and near-future experiments. Markov-Chain Monte Carlo fits to our synthetic 
dataset, which includes a realistic galactic foreground, a plausible model for the signal, and noise consistent 
with 100 hours of integration by an ideal instrument, suggest that a simple four-parameter model that links 
the production rate of Lyman-a, Lyman-continuum, and X-ray photons to the growth rate of dark matter halos 
can be well-constrained (to ^0.1 dex in each dimension) so long as all three spectral features expected to 
occur between 40 < v/MHz <120 are detected. Several important conclusions follow naturally from this 
basic numerical result, namely that measurements of the global 21-cm signal can in principle (i) identify the 
characteristic halo mass threshold for star formation at all redshifts z > 15, (ii) extend z < 4 upper limits on 
the normalization of the X-ray luminosity star-formation rate (L^-SER) relation out to z ~ 20, and (iii) provide 
joint constraints on stellar spectra and the escape fraction of ionizing radiation at z 12. Though our approach 
is general, the importance of a broad-band measurement renders our findings most relevant to the proposed 
Dark Ages Radio Explorer, which will have a clean view of the global 21-cm signal from ^ 40 — 120 MHz 
from its vantage point above the radio-quiet, ionosphere-free lunar far-side. 

Subject headings: early universe — diffuse radiation - epoch of reionization 


1. INTRODUCTION 

The high redshift (z > 6) Universe has become a frontier 
in recent years, as it was the time in which stars, galax¬ 
ies, and black holes first formed, bringing an end to the 
cosmic “dark ages” and initiating the “cosmic dawn.” Pre¬ 
liminary searches for these objects have commenced, pri¬ 
marily with the Hubble Space Telescop e, and have begun 
to find very bright galaxies at z 'y 10 (IZheng et al.l 120121 : 
lEllis et ani2013uOesch et d1120131) . However, it is the in¬ 
dividually faint but overwhelmingly numerous galaxies that 
likely usher in the Epoch of Reionization (EoR), and may 
account for a substantial fraction of the ionizing photons re¬ 
quire d to bring the EoR to a close by z 6 (e.g.- lTrenti et alJ 
120101 IWise etaTI 120141) . If the galaxy lumin osity function 
flatten s considerably at the low-luminosity end (lO’Shea et alJ 
120151) or the formation efficiency of massive Population III 
stars is low, even the James Webb Space Tele scope may strug¬ 
gle to find faint galaxies beyond z ~ 10 (IZackrisson et al.l 
120121) . Measurements of the cosmic microwaye background 
(CMB), in conjunction with independent constraints on the 
ionization and star-formation histories, support a relatiyely 
short EoR, and thus a relatiyely modest galaxy population 
at z ^ 10 (ICollaboration et al.l 120151 iRobertson et al.l IMISl 
iBouwens et al.ll2015l). 

Obseryations at low radio frequencies - corresponding to 
highly redshifted 21-cm “spin-flip” radiation from neutral hy- 
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drogen atoms - are a promising complement to CMB mea¬ 
surements and optical and near-infrared imagin g campaigns 
to constrain the high redshift galaxy population (IMadau et al.l 
119971 IShayer et al.ll 19991 iFurlanetto et al.l l2006^ . Numerous 
efforts to detect the 21-cm background are already underway, 
including both its spatial fluctuations and monopole spectral 
signature. For instance, obseryations with the Precision Ar¬ 
ray for Probing the Epoch of Reionization indicate an X-ray 
heated IGM at z ^ 7.7 — 8.4 (iParsons et al.l 120141 lAli et al.l 
120151 iPober et ^12015h . while the Experiment to Detect the 
Global EoR Signal (EDGES) is so far the only 21-cm exper¬ 
iment to set lowe r limits on the duration of the reionization 
epoch (Az > 0.06: lBowman & Rogersll2010h . 

Targeting the sky-ayeraged (“global”) 21-cm spectrum 
may enable more rapid progress at the highest redshifts 
(z > 10), as it can in principle be detected with a sin¬ 
gle well-calibrated dipole receiyer. Seyeral ground-based 
experiments are underway (e.g.; EDGES, SCTHI, LEDA , 
BIGHORNS: iBowman & RogersI lIoTol iVoytek et al.l Iml 

iGreenhill & Bernardil 120121 ISokolowski et al.ll2015h . though 
space-based obseryatories will be required to probe the cos¬ 
mic daw n at z 30 (e.g. , the Dark Ages Radio Explorer, 
(DARE): lBurns et al.l2012h . as the Earth’s i onosphere reflects 
and refracts radio signals at low frequencies (IVedantham et al.l 
120131 iDatta et ^ 120141) . Interferometers are in principle 
capable of measur ing the global 21-cm signal, so long as 
they are compact (iPresley et al.l 120151 ISingh et aH 120151) or 
can pyercome challeng es in lunar occultation techniques 
(IVedantham et al.ir2014l) . 

No matter the obserying technique, all global experiments 
must cope with the Galactic foreground, which is ^ 10^ — 10® 
times brighter than the cosmological signal in temperature. 
Though strong and sp atially yariable, the Galactic fo reground 
is spectrally smooth (Ide Oliyeira-Costa et al.ll2008h . in con¬ 
trast to the expected high redshift signal which is spatially 
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invariant but spectrally complex. It is this spectral structure 
that should enable one to distinguish foreground from signal 

(IShaver et alJll99^lGnedin & Shaverll2004t lFurlanett3l2006l : 

iPritchard & Loebll2010h . especi ally If one observes multiple 
(semi - ) independent sky reg ions (iHarker et al.l2012tlLiu et aP 
12013l ISwitzer & Liull2()T^ . Though current constraints do 
not rule out unresolved spectral structure at the level of the 
high redshift signal, theoretical argu ments favor a smooth 
foreground at the re levant frequencies (iPetrovic & Ohll20Tll : 
iBernardi et al.ll20f5ll . Extragalactic point sources are another 
foreground, but for experiments with broad beams, their com¬ 
bined contribution averages into anot her diffuse spectrally 
smooth foreground (IShaver et al.lll999h . 

Given that the global 21-cm signal is an indirect probe of 
high-z galaxies, some modeling is required to convert ob¬ 
servational quantities to constraints on the properties of the 
Universe’s first galaxies. Though numerous studies have per¬ 
formed forwa rd modeling to predict the strength of the g lobal 
21-cm signal (iFurlanettolUOCibt [Pritchard & Loebll201()l) . few 
have attempted to infer physical parameters of interest from 
synthetic datasets. Such forecasting exercises are incredibly 
useful tools for designing instruments and planning observ¬ 
ing strategies, as they illuminate the mapping between con¬ 
straints on observable quantities and model parameters of in¬ 
terest. Both Fisher matrix and Markov Chain Monte Carlo 
(MCMC) approaches h ave been employed by the power spec¬ 
trum community (e.g., iPober et ani2014t iGreig & MesingeJ 
12015h . the latter proyiding a powerful generalization that does 
not require the assumption of Gaussian errors or perfect re- 
coyery of the maximum likelihood point. 

Most work to date has instead focused on forecasting con¬ 
straints on phenomenological para meters of intere st, e.g., the 
timing and duration of reionization (iFiu et al.l2013[) . the depth 
and width of the de ep minimum expecte d near ^70 MHz 
prior to reionization (IBernardi et al.lUOT^ . or all three spec- 
I ral features predicted to occur between 40 5. V /MHz <120 
(IPritchard & FoebI 1201 Ot IHarker et al.l 1201^ . These spec¬ 
tral “turning points” in the global 21-cm spectrum proyide a 
natural basis for parameter forecas ting, as they persist oye r 
large ranges of parameter space (IPritchard & FoebI l2010h . 
can be extracted from the foreground with realistic instru¬ 
ments and integration time s (at least under the assumption 
of a negligible ionosphere ; IHarker et al.l 120121: iPresley et al.l 
120151: IBernardi et al.l iMIsI) . and can be extracted from the 
foreground eyen when the ir positions are not used as the pa¬ 
rameters of a signal model (IHarker et al.l2015. in prep.ll . They 
can also be interpreted fairly robustly in terms of the physical 
properties of the IGM , at least in simple two-phase models 
(iMirocha et al.ll20f3l hereafter Paper I). Giyen the yiability of 
the turning points as “products” of global 21-cm signal ex¬ 
traction pipelines, we will use them as a launching point in 
this paper from which to explore the prospects for constrain¬ 
ing astrophysical parameters of interest with obseryations of 
the global 21-cm signal. Importantly, we will consider all 
three turning points simultaneously, rendering our findings 
particularly applicable to DARE, whose band extends from 
40 < v/MHz < 120 in order to maximize the likelihood of 
detecting all three features at once. 

This paper is organized as follows. In Section 2 we outline 
our methods for modeling the global 21-cm signal and pa¬ 
rameter estimation. Section 3 contains our main results, with 
a discussion to follow in Section 4. In Section 5, we sum¬ 
marize our results. We use the most up-to-date cosmological 
parameters from Planck throughout (last column in Table 4 of 


ICollaboration et al1l2015l) . 

2. NUMERICAL METHODS 

In order to forecast constraints on the properties of the first 
galaxies, we will need (1) a model for the global 21-cm signal, 
(2) estimates for the precision with which this signal can be 
extracted from the foregrounds, and (3) an algorithm capable 
of efficiently exploring a multi-dimensional parameter space 
to find the best-fit model parameters and their uncertainties. 
The next three sub-sections are deyoted to describing these 
three pieces of our pipeline in turn. 


2.1. Physical Model for the Global 21-cm Signal 

Our approach to modeling the global 21-cm signal is sim- 
ilar to that presented in seyeral other pub l ished works (e.g., 
Barkana & FoebI 120051 iFurlanettol 120061: IPritchard & FoebI 
20101: lMirochall20l4 . so we will only discuss it here briefly. 
The primary assumption of our model is that the radiation 
backgrounds probed by the turning points are generated by 
stars and their byproducts, which form at a rate proportional 
to the rate of baryonic collapse into dark matter haloes. That 
is, we model the star-formation rate density (SFRD) as 

= (1) 

where f^ou = /coiK^min) is the fraction of matter in collapsed 
halos with yirial temperatures greater than Tmin, pj^ is the 
mean baryon density today, and /* is the star formation ef¬ 
ficiency. We use a fixed Tmin rather than a fixed because 
it proyides physical insight into the processes goyerning star- 
formation, as one can easily identify the atomic and molec¬ 
ular cooling thresholds of ^ 10^ and ~ 500 K. Note that a 
fixed yalue of Timn results in a time-dependent mass thresh¬ 
old, M^nin- 

In order to generate a model realization of the global 21-cm 
signal, we must conyert star-formation to photon production. 
Giyen that the three spectral turning points probe the history 
of ionization, heating, and Fy-a emission, we will split the 
production of radiation into three separate bands: (1) from 
the Fy-a resonance to the Fyman-limit, 10.2</rv/eV < 13.6, 
which we refer to as the Fyman-Werner (FW) band despite 
its inclusion of photons below 11.2 eV, (2) hydrogen-ionizing 
photons, with energies 13.6 < /rv/eV < 24.4, and (3) X-rays, 
with energies exceeding 0.1 keV. Fach radiation background 
is linked to the SFRD through a scaling parameter which 
conyerts a rate of star formation to a rate of photon (or energy) 
production). The Fy-a background intensity is then giyen by 


, . _ c (l-bz)^ ^LW-Q <i/coll 

“ ^ 47t H{z) AVa^'’ dt ’ 

where AVa = Vll — Va, and the ionization rate by 


rH,= 


S ^0 ^/coll 

dt 


( 2 ) 

(3) 


The rate of X-ray heating is defined instead in terms of an 
energy per unit star-formation, i.e.. 


= fxMCX^xPb^^^ (4) 

where cx is the normalization of the L^-SFR relation, which 
we take to be cx = 3.4 x lO^*’ erg s^* {Mq yr)^' following 
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iFurlanettol (l2006iri Note that we have absorbed /* into the ^ 
parameters, i.e. 

^LW = Nlw/* (5) 

^ion = Nion/*/esc (6) 

'^x=f*fx, (7) 

where Nlw and Nion are the number of LW and ionizing pho¬ 
tons emitted per stellar baryon, /esc is the escape fraction of 
ionizing radiation, and fx scales the Lx-SFR relation. 

Given ex, Fhi, and the Fy-a background intensity, Ja, we 
can evolve the ionization and thermal state of intergalactic gas 
in time, and compute the sky-averaged 21-cm signal, i.e., the 
differential brightness temperature of HI relative to the CMB, 
via (e.g., lFurlanett^l2006lf 1 


57;-27(1-1,) 


V 0.023 J 


f 0.15 i+^ v^7i 
\a^,oh^ 10 J V Ts/’ 

( 8 ) 


where Ty is the CMB temperature, x, is the volume-averaged 
ionization fraction. 


Figure [T] also shows that will be difficult to constrain 
using global signal data at these frequencies, as even factor of 
10 changes lead only to small changes in the signal (at v > 100 
MHz), whereas factor of 10 changes in ^x and ^lw are — 50 
mK effects. There are also clear degeneracies between 
and the ^ parameters. Exploring those degeneracies and de¬ 
termining the prospects for constraining each parameter inde¬ 
pendently are our primary goals in this work. The results will 
in large part depend on how accurately the signal can be re¬ 
covered from the foregrounds, which we discuss in the next 
subsection. 

Before moving on to signal recovery, it is worth reiterating 
that our approach cannot be used to constrain the normaliza¬ 
tion of the SFRD, since we have absorbed the star-formation 
efficiency into the ^ parameters. However, we can constrain 
the rate-of-change in the SFRD, as it is uniquely determined 
by T^nin. It is illustrative to quantify this using an effective 
power-law index 


aeff(z) 


^logPH.(z) 
d\og{\ -l-z) ’ 


( 11 ) 


Xi = Q\\u + {^ — Qm\)xe (9) 


with Qmi representing the volume-filling factor of HII re¬ 
gions and Xg the ionized fraction in the bulk IGM. Tj is the 
excitation or “spin” temperature of neutral hydrogen, which 
quantifies the number of hydrogen atoms in the hyperfine 
triplet and singlet states. 


J Ty ' -\-XcTfy ' -fxaT),, ' 

1 + Xt + Xfx 


( 10 ) 


where Tx is the gas kinetic temperature, Ta — 7/ (lFieldll958l) . 
and h and the H’s take on their usual cosmological meaning. 
We compute the collisional coupling coe fficient, Xn, by in- 
terpolating between the tabulated values in IZveelmanl (l2005h 
with a cubic spline, and takexa = 1.81 x 10**7a/(l -|-z). We 
perform these calculations using the Accelerated Reionization 
Era Simulations (ARES) codfl which is the union o f a 1- 
D radiative transfer code developed in iMirocha et al.l (12012h 
and un iform radi ation backgroun d code described in iMirochal 
dSoIl. See §2 of iMirochal (EoTl) for a more detailed descrip¬ 
tion of the global 21-cm signal modeling procedure. 

Figure [T] shows our reference model (properties of which 
are listed in Table [D. and the modulations in its structure that 
occur when varying Ti^iin, ^x, ^ion, and ^lw- It is immediately 
clear that 7/iin affects the locations of all three turning points, 
whereas each ^ parameter affects at most two. We should 
therefore expect that in principle, an observation containing 
all three features will have the best chance to constrain Tmin, 
though this will be complicated at the lowest redshifts where 
dfcow/dt becomes a weaker function of T^m (see bottom panel 
of Figure E). 


^ IFurlanettol ZOOfl) comput ed thi s value by extrapolating the 2-10 keV Lx- 
SFR relation of IGrimm et alj <200311 above 0.2 keV, assuming an unabsorbed 
a = 1.5 power-law spectrum. Our reference value of 4x = ftfx = 0.02 
(which is dimensionless, unlike ^lw and ^ion) is chosen to match recent anal- 
vses in the 0.5-8 ke V band, which find cx = 2.6 x 10^^ erg s“' {Mq yr)“* 
<Mineo et al.l20123l 

* Though this expression does not explicitly depend on the density of hy¬ 
drogen gas (it assumes gas is at the cosmic mean density, 8 = 0), the global 
21-cm signal may still be sensitive to the density provided that fluctuations 
in the ionization and/or spin temperature of FlI gas are correlated with its 
density. 

® https;//bitbucket.org/mirochaj/ares;vO.1 


which enables a straightforward comparison with empirical 
models, which are often power-laws, i.e., p*(z) (1 in 

which case a = tteff = constant. The tteff (z) values of our /coii 
model are independent of /* so long as t //*/dt = 0, and gener¬ 
ally fall within the (broad) range of values permitted by obser¬ 
vation s of high-z galaxies (lOesch et al.ll2013l : lRobertson et alJ 
120151) . as shown in the top panel of Figure |2] 

2.2. Signal Extraction 

In order to fit a physical model to the turning points of the 
global 21-cm signal, we require best-fit values for the turn¬ 
ing point positions and estimates for uncerta intie s. To do this, 
we build on the work of lHarker et alJ (1201 2^ and iHarker et alJ 
dloK in prepJ) . who introduced a MCMC technique for fit¬ 
ting global 21-cm signal data. The basic approach is to simul¬ 
taneously fit a model for the galactic foreground, the global 
21-cm signal, and in general, parameters of the instrument 
(e.g., its response as a function of frequency), assuming some 
amount of integration time, hnt, and the number of indepen¬ 
dent sky regions observed, iVsky. The foreground is modeled 
as a polynomial in logv — logT space, while the astrophy si- 
cal signal is modeled as either a spline (iHarker et al.ll20121) or 
series of tanh functions, 

A- f 

A(z) = ^{1+tanh[(zo-z)/Az]} , (12) 

that represen t 7a (z), 7k(z), and x/z) (IHarker et alJ 
l20K in prepJ) . The free parameters of the tanh model 
are the “step height,” Aj-gf, pivot redshift, zo, and a duration, 
Az. 

The tanh approach to modeling the global 21-cm signal 
was chosen for numerous reasons. First and foremost, it was 
chosen as a computationally efficient substitute for more ex¬ 
pensive, but physically-motivated models like those investi¬ 
gated in this paper. Some alternativ e intermediaries include 
the ‘turning points ’ parameterization (iPritchard & Foebl2010l : 
IHarker et alJl2012h or models that tre at the absorption feature 
as a Gaussian (iBernardi et al.ll20L5l) . Both are comparably 
cheap computationally, but cannot capture the detailed shape 
of physical models. Perhaps most importantly, the spline and 
Gaussian models are purely phenomenological, making them 
difficult to interpret in terms of IGM or galaxy properties and 
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Figure 1. Illustration of the basic dependencies of the global 21-cm signal. The black line is the same in each panel, representing our reference model (see Tables 
[T]and|2), while all solid green (blue) lines correspond to a factor of 2 increase (decrease) in the parameter noted in the u pper left comer, and dash ed lines are 
factor of 10 changes above and below the reference value. The right half of the figure is qualitatively similai* to Figure 2 of IPritchard & Lo^ 420101) . though our 
reference values for the ^ parameters are different, as are our cosmological parameters, leading to quantitative differences. The dotted lines show the maximum 
allowed amplitude of the signal (i.e., the amplitude of the signal when Xj — 0 and Ts >> Ty), and the minimum allowed amplitude of the signal (set by assuming 
Ts — Tk — Tad(z), where Tad is the gas temperature in an adiabatically-cooling Universe). Because we refer to the spectral features as turning points B, C, and D 
throughout the paper, we annotate them in the lower left panel for reference. 


Table 1 

Reference Model Properties and Simulated Constraints 


Quantity 


Reference Value 

EMI 

EM2 

Description 

Vb [MHz 


47.4 

46.99 ±0.74 

47.08 ±0.60 

Onset of Ly-a coupling 

Vc [MHz 


71.0 

70.95 ±0.20 

70.96 ±0.15 

Onset of heating 

Vd [MHz 


111.4 

110.9±5.0 

109.2±3.5 

Beginning of reionization 

sri(vB) 

mK 

-4.4 

n/a 

n/a 

Depth when Ly-a coupling begins 

87i,(vc) 

mK 

-124.8 

-122.6±5.0 

-121.7±3.7 

Depth of absorption trough 

S7i(VD) 

mK 

19.2 

17.20±4.5 

19.88±1.7 

Height of emission feature 

2rei 


535 

n/a 

n/a 

Midpoint of reiomzation 

Xe 


0.066 

n/a 

n/a 

CMB optical depth 


Note. — Observational properties of our reference model (solid black lines in FigureQ, and the best-fit and uncertain¬ 
ties for each extraction model (EM) we consider. Subscripts indicate different turning points, i.e., the cosmic dawn feature 
when the Wouthuysen-Field effect first drives Ts to Tk (B), the absorption trough, which indicates the onset of heating 
(C), and the beginning of reionization (D). EMI and EM2 differ in the number of independent sky regions assumed (1 vs. 
2), and in the complexity of the foreground model (3rd vs. 4th order polynomial), which leads primarily to a more robust 
detection of turning point D for EM2. All errors are 1 — rt, and correspond directly to the diagonal elements of the turning 
point covariance matrix. 


thus incapable of incorporating independent prior information 
on e.g., the ionization or thermal history. The tanh approach, 
on the other hand, can mimic the shape of typical global 21- 
cm signal models extremely well, and can be immediately re- 
la ted to physical properties of t he IGM. 

iHarker et al.l (12015. in r)rep.h presented a suite of calcula¬ 
tions spanning the 2-D parameter space defined by Vsky = 
{1,2,4,8} and tint = (100,1000}. In the tint = 1000 hr cal¬ 
culations, confidence contours of the turning point positions 


narrowed enough to reveal subtle biases in their recovered po¬ 
sitions, which led to biases in constraints on physical proper¬ 
ties of the IGM as well. These shifts were interpreted to be 
due to degeneracies between the signal and the foreground 
at high frequencies, as they could be mitigated by using a 
more sophisticated foreground model or, at the expense of 
losing information from turning point D, simply by truncat¬ 
ing the bandpass at 100 MHz. However, even with unbiased 
constraints on the turning point positions, biases in the IGM 
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Figure 2. Bottom: Rate of collapse onto halos above a given virial temper¬ 
ature thres hold , Tmin, scaled to the Hubble time. Top: Effective power-law 
index (Eq. E) a s a function of reds hift for each t//coii /dt model. Empirical 
power-laws from lOesch et all 120121) are ove rlaid for comparison, a s well as 
the best-fit 4-parameter SFRD model used in IRobertson et all 120151) . 

properties persisted, likely because the tanh is not a perfect 
match in shape to the physical model injected into the syn¬ 
thetic dataset. Though we suspect that biases in constraints 
on the turning points are a persistent feature in datasets em¬ 
ploying the tanh, turning point constrain ts from the fjnt = 100 
hr cases in (iHarker et al.ll2015. in prep.l) are broad enough to 
hide such biases. In using the fim = 100 hr results, we should 
then expect to be able to obtain unbiased constraints on the pa¬ 
rameters of our physical model in the present work. We also 
only analyze cases using one and two sky regions, for which 
model s election will b e more immediately tractable computa¬ 
tionally (lHarkeill2015h . 

Now, back to the simplest model of IHarker et alJ 
(12015. in prep.h (EMI in Table [T]i. This calculation assumed 
a single sky region, 100 hours of integration, and a third- 
order logv — logT polynom ial for the galactic foregrouncP^. 
IHarker et alJ (12015. in prep.h investigated the generic case 
of an idealized instrument (a flat 85% response function), 
though this could easily be modified to enable forecast¬ 
ing for non-ideal instruments. The foreground and astro- 
physical signal were simultaneously fit using the parallel- 
tempering sampler in the pub licly available EMCEE codJ^ 
(lEoreman-Mackev et al.ll2M3ll . a PYTHO N implementation of 
the af fine-invariant MCMC sampler of iGoodman & Wear^ 
(1201 Oh . from which constraints on the positions of the turning 
points followed straightforwardly. The errors on the turning 
points are in general not Gaussian, and are often correlated 
with one another, though for the purposes of our fitting, we 
approximate the errors as 1-D independent Gaussians since 

It seems lilcely that in practice a higher order polyn omial will be needed 
to fit out instrumental effects te.g.. lBernardi et M2015h . However, here our 
synthetic datasets contain foregrounds with no structure beyond a polynomial 
of order three (or four), meaning third and fourth order polynomials can fit 
the foreground perfectly (by construction). 

http://dan.iel.fm/emcee/current/ 
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Figures. Comparison of EMI and EM2 for turning point D, the point at 
which they differ most substantially (see Table [T). In blue and red, solid 
(dashed) contours denote 68% (95%) confidence regions for EMI and EM2, 
respectively. The dotted black line shows the saturated limit, in which J, = 
0 and Ts ~ Tk ^ Ty, while the dashed vertical line denotes the position of 
turning point D in our reference model (solid black curve). Note that the 
EMI error ellipse for turning poin t D extends to ^ 130 MHz, beyond the 
edge of the bandpass considered in lHarker et all <2015. in prep.l) , though the 
2 — c upper limit for EM2 is within the assumed band, at Vd ~ 117 MHz. 


covariances are likely to depend on the choice of signal pa¬ 
rameterization. In addition to EMI, we also investigate the 
results of a fit using a 4* order log-polynomial for the fore¬ 
ground model, which we refer to as EM2. Table [T] summa¬ 
rizes the different signal extraction models, while Eigure [3] 
illustrates the primary difference between the two extraction 
models graphically. 

2.3. Parameter Estimation 

With a physical model for the global 21-cm signal (|2Jl 
and a set of constraints on the turning point positions ( ^2.2b . 
we then explore the posterior probability distribution function 
(PDE) for the model parameters, 0, given the data, (D. That is, 
we evaluate Bayes’ theorem, 

p(e|(D)oc i:((D|e)(p(e). (13) 

The log-likelihood is given by 

logZ(D|e)~-5:M5^ (14) 

where p, is the “measurement” with errors a, (i.e., the values 
listed in columns 3 and 4 of Table [Til, and x(0) represents 
a vector of turning point positions extracted from the model 
global 21-cm signal generated with parameters 0. This “two- 
stage approach” to fitting the gl obal 21-cm s i gnal - the first 
stage having been conducted bv iHarker et al.l (12015. in prep.1) 
- is much more tractable computationally than a direct “one- 
stage” fit to a mock dataset using a physical model. Note 
that the brightness temperature of turning point B, 87), (vb), 
is tightly coupled to its frequency, so we are effectively only 
using 5 independent data points in our fits. 

To explore this four-dimensional space, we use EMCEE. We 
assume broad, uninformative priors on all parameters (listed 
in Table |2]i, but note that our physical model implicitly im¬ 
poses three additional constraints on the astrophysical signal: 
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Table 2 

Parameter Space Explored 


Parameter 

Description 

Input 

Min 

Max 

Tm,„ (K) 

Min. virial temp, of star-forming haloes 

10'* 

100 

lO^-’ 

^LW 

Ly-a efficiency 

969 

10 

lOi^ 


X-ray efficiency 

0.02 

10--* 

10** 

^io„ 

Ionizing efficiency 

40 

10-* 

10^ 


Note. — Parameter space explored for results presented in 0 The first two columns 
indicate the parameter name and a brief description, the third column is the “true value” of 
the parameter in our reference model, while the last two columns indicate the bounds of the 
priors for each parameter, all of which are assumed to be uninformative, i.e., modeled as 
uniform distributions between the minimum and maximum allowed values. 


1. We neglect exotic heat sources at high-z, which con¬ 
fines turning point B to a narrow “track” at v < 50 MHz. 

2. We assume that the Universe cannot cool faster than 
the Hubble expansion, which sets a redshift-dependent 
lower limit on the strength of the absorption signal 
(lower dotted curve in all panels of Figure[Tli. 

3. We assume the mean density of the IGM we observe is 
the universal mean density, i.e., it has 8 = 0, which pre¬ 
vents the signal from exceeding the “saturated limit,” 
in which Tg ^ T^and x, = 0 (upper dotted curve in all 
panels of Figure [ill. 

Our code could be generalized to accommodate exotic heating 
models, though this is beyond the scope of this paper. Bul¬ 
lets 2 and 3 above are manifestly true for gas at the cosmic 
mean density (via Equation [^, though imaging campaigns 
will likely see patches of IGM whose brightness tempera¬ 
tures exceed (in absolute value) these limits, owing to over¬ 
densities 5 > 0. 

For all calculations presented in this work, we use 384 
walkers, each of which take a 150 step burn-in, at which point 
they are re-initialized in a tight ball centered on the region of 
highest likelihood identified during the burn-in. We then run 
for 150 steps more (per walker), resulting in MCMC chains 
with 57,600 links. The mean acceptance fraction, i.e., the 
number of proposed steps that are actually taken during our 
MCMC runs, is ^ 0.3. The runs are well-converged, as we see 
no qualitative differences in the posterior distributions when 
we compare the last two 10,000 element subsets of the full 
chain. 

3. RESULTS 

Each MCMC fit yields 57,600 samples of the poste¬ 
rior distribution, which is a 4-dimensional distribution in 
{7mm,^LW,^x,^ion} Space. However, we also analyze each 
realization of the global 21-cm signal on-the-fly as the 
MCMC runs, saving IGM quantities of interest every Az = 0.1 
between 5 < z < 35, as well as at the turning points. To build 
upon the analytical arguments presented in Paper I, which 
provided a basis for interpreting the turning points in terms 
of IGM properties, we start with an ana lysis of the inferred 
IGM properties at the turning points in ^3.11 deferring a full 
analysis of the IGM history to future work. Readers inter¬ 
ested only in the co nstra ints on our four-parameter model can 
proceed directly to ^3.21 

3.1. Constraints on the Intergalactic Medium 


We begin by showing our mock constraints on properties 
of the IGM at the redshifts of turning points B, C, and D in 
EigureslHIS] and [6] respective! j{3- 
Because turning point B primarily probes the Ly-a back¬ 
ground, we focus only on its ability to constrain Ja in Eigure 
m The input value is recovered to 1—0, with relatively tight 
error-bars limiting Ja to within a factor of 2. Use of EM2 has 
little effect on this constraint as its main improvement over 
EMI is at frequencies v > 100 MHz. 



Figure 4. Constraints on the background Ly-a intensity at the redshift of 
turning point B, in units of Jii — 10“^*erg s“^ cm“^ Hz“* sr“^. Dotted 
vertical line shows the input value, which occurs at z = 29 in our reference 
model. The black histogram is the constraint obtained if using EMI, while 
the analogous constraint for EM2 is shown in green. 

Eigure [5] shows constraints on the Ly-a background and 
thermal history at the redshift of turning point C. In the ~ 90 
Myr separating turning points B and C, the Ly-a background 
intensity, Ja, has risen by a factor of ~ 350, though we still 
constrain its value to within a factor of ~ 2 (panel f). The 
IGM temperature is limited to 9 < 7k/K< 11, and would 
otherwise be ~ 7.4 K at this redshift in the absence of heat 
sources. There are noticeable degeneracies in the 2-D PDEs, 
which are not necessarily obvious intuitively. 

Let us first focus on the anti-correlations in the Ja-£x and 
Ja — Tk planes (panels a and d in Eigure |5]l. Eor this exercise 
- and those that follow - it will be useful to consider slight 
excursions away from our reference model. We can see from 
the lower right panel of Eigure [T] that a small increase in ^lw 
will shift turning point B to slightly higher redshifts (lower 
frequencies) holding all other parameters fixed. Turning point 
C will also occur earlier than in our reference model (since a 
stronger Ly-a background can couple 7$ to 7k more rapidly) 
and be deeper, since there has been less time for X-rays to heat 
the IGM, leading to increased contrast between the IGM and 
the CMB. Panels (a) and (d) in Pigure|5]now make sense: the 
anti-correlations in the 7a-ex and Ja — Tk planes arise because 
measurement errors permit slight excursions away from the 

Note that our choice to derive constraints on the IGM parameters at the 
turning points, rather than at a fixed series of redshifts, is in part responsible 
for some of the behavior in Figures and [6] This seemed a natural choice 
given that we fit our model to the turning points alone. Constraints at an 
arbitrary redshift could be extracted in future studies, for example to compare 
to independent measurements that do not coincide with the redshifts of the 
turning points. 
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reference model, which if achieved through enhancements to 
^LW, shift turning points B and C to slightly earlier - and thus 
cooler - times. 

One could also counteract a mild increase in ^lw with a cor¬ 
responding increase in which enhances heating and thus 
leads to shallower absorption troughs. However, increasing 
^x shifts turning point C to shallower depths and lower fre¬ 
quencies, thus exacerbating the leftward shift caused by larger 
values of ^lw- As a result, Emin would also need to be in¬ 
creased in order to delay the onset of Wouthuysen-Field cou¬ 
pling and heating. Indeed, we will find this series of positive 
corr elati ons among the physical parameters of our model in 
the ^ 

Before moving on to the IGM constraints associated with 
turning point D, we note that the correlation between Tx and 
ex (panel b) is simply because Tk°^ f exdz, and ex is mono¬ 
tonic. Also, apparently the improvement at the highest fre¬ 
quencies offered by EM2 also acts to slightly bias constraints 
on Ja and ex relative to their input values. Referring back to 
Figure |3l we do see a slight bias in the EM2 PDF for turning 
point D toward larger amplitude, which would require more 
rapid heating at earlier times. In fact, this is precisely the 
sense of the bias we see in Figure |5] slightly larger values of 
ex at turning point C, and a corresponding downward shift in 
Ja as described above. 



Figure 5. Constraints on the Ly-ct background intensity, IGM temperature, 
heating rate density at the redshift of turning point C. The heating rate den¬ 
sity, e^, is expressed in units of erg s“* cMpc“^, while Ja is once again 
expressed in units of 721 = 10“^*erg s“^ cm“^ sr“*. Dotted vertical 

lines show the input values, which occur at z— 19 in the reference model. 
Open contours are those obtained with EMI (68% and 95% confidence re¬ 
gions in solid and dashed curves, respectively), while filled contours are the 
constraints obtained by EM2 (68% and 95% confidence regions in blue and 
green, respectively). The color-scheme along the diagonal is the same as in 
Figure[4] with EMI (EM2) curves in black (green). 


And finally. Figure |6] shows constraints on the ionization 
and thermal histories at the redshift of turning point D, which 
occurs at z = 11.75 in our input model. The behavior here 
is complex, as the signal is not yet saturated (i.e., Ts ^ Ty is 
a poor approximation) and the mean ionized fraction is non¬ 
zero (i.e., Qmi ^ 0.2). This means the global 21-cm signal 
depends on both the ionization history and the thermal his¬ 
tory, which we may parameterize in terms of the volume fill¬ 
ing factor of ionized gas, Qmi, the IGM temperature, Tk, and 


their time-derivative^ Fhi and ex- We may, however, ne¬ 
glect the Fy-a history at this stage, since Ts ~ Tk is accurate 
to high precision, rendering any constraints on Ja completely 
parameterization-dependent (i.e., Ja can be anything, so long 
as it is large enough to drive Ts Tk). 

It is once again useful to consider excursions away from the 
reference model. At fixed thermal history, a slight increase in 
^ion will act to decrease the amplitude of turning point D and 
shift it to slightly higher redshift. With less time to heat the 
IGM between turning points C and D, the IGM is cooler at the 
redshift of turning point D in this scenario and as a result, the 
emission signal is weaker than that of our reference model. 
This line of reasoning explains the anti-correlations between 
the ionization and thermal history parameters in Figure|6l As 
in Figure |5] positive correlations occur by construction, since 
state quantities like Qmi and Tk are just integrals of Fhi and 
ex, which are both monotonically increasing with decreasing 
redshift. 

The advantages of EM2 over EMI are also clear in Figure|6l 
This improvement occurs because EMI does not detect turn¬ 
ing point D with significance away from the saturated limit or 
within the assumed band (v < 120 MHz), whereas the EM2 
fit does both at the > 2 —a level. Perhaps most notably, this 
leads to a strong detection of the early stages of reionization 
(0.12 < Qmi < 0.29 at 2 — a; green PDF in panel j). 

Fastly, we note that although the amplitude of the signal 
is set by x;, a volume-averaged ionized fraction, we only 
show constraints on Qmi, as Xg never reaches values above 
^ 10^^ at z > 10 in any of our calculations. As a result, 
it has a negligible impact on the ionization history. How¬ 
ever, even mild ionization of the bulk IGM enhances the ef¬ 
ficiency of heating rather substantially since the fraction of 
photo-electron energy deposited as heat (as opposed to ion¬ 
ization or excitation) is a strong function of the electron den¬ 
sity ( e.g. IShull & van Steenberglll985l: iFurlanetto & Stoe\^ 
l201(]h . which means the value of Xg can have a considerable 
effect on the thermal history of the IGM. Our choice of a mean 
X-ray photon energy of /zVx =0.5 keV, in lieu of a detailed 
solution to the radiative transfer equation, drives t his result. 
More detailed calculations that solve the RTF (e.g. iMirochal 
12014) could enable scenarios in which the bulk IGM is ion¬ 
ized substantially prior to the overlap phase of reionization, 
which could have interesting observational signatures. We de¬ 
fer a more detailed treatment of this effect to future work. 

3.2. Constraints on the Physical Model 

Our main results are illustrated in Figures|7][8] and|9] which 
analyze the full 4-D constraints on our reference model and 
the implications for UV and X-ray sources, respectively. In 
this section, we’ll examine each in turn. 

It is perhaps most intuitive to begin with the panels along 
the diagonal of Figure |7] which show the marginalized 1-D 
constraints on the parameters of our reference model. As 
predicted, given its broad-band influence on the signal, Tmin 
(panel d) is most tightly constrained, with la error bars of or¬ 
der ^ 0.05 dex. Therefore, an idealized instrument observing 
a single sky region for 100 hours can rule out star-formation 

Although we use the symbol F, we caution that our values should not 
be compared to extrapolations of constraints on Fm from the Ly-a forest at 
2 < 6. The latter is a probe of the meta-galactic ionizing background (i.e., 
large-scale backgrounds), whereas our values of F probe the growth rate of 
ionized regions, and thus should be considered a probe of radiation helds near 
galaxies. A more detailed cosmological radiative transfer treatment could in 
principle reconcile the two tracers of ionizing sources. 
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Qhii logio^x 


Figure 6. Constraints on the volume-filling factor of ionized gas, volume-averaged ionization rate, IGM temperature, and heating rate density at the redshift of 
turning point D, which occurs at z = 11.75 in our reference model. Open contours ai‘e 68% (solid) and 95% (dashed) confidence regions for EMI, while filled 
contours show the results from EM2, with 68% and 95% confidence regions shown in blue and green, respectively. Input values are denoted by black dotted lines 
in each panel. 


in molecular halos (onto which gas collapses more slowly; see 
Figure|2]i, at least at levels sufficient to affect all three turning 
points. Errors on and ^lw are comparable (panels i and 
j), though the positive error-bars are larger at ^ 0.1 dex. The 
errors on ^ion are more asymmetric, at -fO. 1/ — 0.2 dex (panel 
g)- 

Strong degeneracies are also apparent, particularly in pan¬ 
els (a), (b), and (h), which show 2-D constraints in the linin- 
^LW, 7min-^x, and ^LW-^x planes, respectively. The first two 
are straightforward to understand. An increase in ^lw drives 
an enhancement in Ly-a production per unit star-formation, 
which can be counteracted by a reduction in the star-formation 
rate density. In our modeling framework, a reduction in the 
SFRD is achieved by increasing Tmin, confining star formation 
to more massive and thus more rare halos. If /* were allowed 
to vary, it too could limit the SFRD, though the change would 
be systematic, whereas varying Tmin affects both the normal¬ 
ization and the redshift evolution. The same line of reasoning 
explains the relationship between Tmin and 

The ^LW-^x degeneracy is slightly more complex. An in¬ 
crease in ^Lw seeds a stronger Fy-a background (holding Tmin 
fixed), which in turn shifts turning point B to lower frequen¬ 
cies (see the lower right panel of Figure [TJ, which measure¬ 


ment error permits to some degree. This will result in a deeper 
(and earlier) absorption trough unless is increased, causing 
a shallower trough (see upper right panel of Figure [U. Once 
again, measurement error sets the degree to which an increase 
in tx c an compensate for an increase in ^lw- As discussed 
in q3.1l slight excursions in ^x cannot completely correct for 
changes in ^lw. and will require changes in T^m as well, es¬ 
pecially if the measurement errors are small. In the limit of 
very large error-bars, however, confidence contours would not 
close and instead we would have large “bands” through pa¬ 
rameter space, signifying an insurmountable degeneracy be¬ 
tween two parameters. Our results indicate that observations 
of a single sky region for 100 hours, albeit with an idealized 
instrument, are precise enough to close these contours, and 
recover all input values to within 1 — o confidence. We will 
revisit this claim in ® 

At this stage it may seem like we have just traded con¬ 
straints on one set of phenomenological parameters (the pa¬ 
rameters of the tanh model; Equation fT2li for another (r^in 
and the ^’s). However, if we assume that ^lw and ^ion probe 
the same stellar population, their ratio is independent of the 
star-formation history (which is set by /* and T^m, which we 
assume are time-independent), and thus constrains the spec- 
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Figure 7. Constraints on our 4-parameter reference model. Filled contours in the interior panels are 2-D marginalized posterior PDFs with 68% confidence 
intervals shaded blue and 95% confidence regions in green. Panels along the diagonal are 1-D marginalized posterior PDFs for each input parameter, with 1-G 
asymmetric error-bars quoted, as computed via the marginalized cumulative distribution functions. Dotted lines denote the input values of our reference model 
(Table^. Bins of width 0.05 dex are used in each panel. Annotated best-fit values and eiTor bars along the diagonal ai‘e those associated with EM2. 


tral energy distribution (SED) of galaxies modulo a factor of 
the escape fractioiQ i.e. (following Equations[5]and[6l), 


^LW 


Vi. 


Vlw 


/esc 


(15) 


To compute Vion/A/w, we take model s pectra l energy dis¬ 
tributions directly from iLeitherer et al.l (Il999h . We focus 
on those assuming an instantaneous burst of star-formation 
with nebular emission included (their Figures 1, 3, and 5), 
and find the cumulative number of photons emitted in the 
LW and hydrogen-ionizing bands, which typically plateaus 
around ^ 20 Myr after the initial burst. The results, as a func¬ 
tion of metallicity and stellar initial mass function (IMF), are 
shown in the left panel of Figure 0 While the values of Vion 
and Vlw vary by factors of ~ 2 over the metallicity range 
0.001 < Z/Zq < 0.04, their ratio changes by only ^ 0.05 over 
this same interval in metallicity. The more important depen¬ 
dence is on the stellar IMF: a standard Salpeter IMF, with 
OCimf = 2.35 and an upper mass cutoff of Mcut = 100 Mq, 
yields 0.25 < Vion/Vuw < 0.3 for all lO^^ < Z/Zq < 0.04 


We assume the escape fraction of LW photons is 100%, though in reality 
this is only likely to be true in the smallest halos (e.g. JKitavama et alJ2004l) . 
For simplicity, we neglect this complication and defer a more thorough treat¬ 
ment to future work. 


(blue circles in the left panel of Figure^, whereas mass func¬ 
tions with fewer massive stars, whether that be achieved with 
steeper power-law indices (ttiMF = 3.3; green squares in Fig¬ 
ure or by reducing the upper cutoff (M.ut = 30 Mq; red 
triangles in Figure[^, yield 0.07 < Von/A/w ^ 0.12. 

In the right panel of Figure[8] we compare our constraints in 
the ^LW — ^ion plane with the Z = 0.04 Zq stellar population 
models described above. The red, green, and blue bands in 
the right panel correspond to the stellar population models 
denoted by filled points of the same color in the left panel. We 
also show the case of a pure 50,000 K blackbody spectrum in 
the cross-hatched region. The width of each band corresponds 
to a factor of two change in the escape fraction, 0.1 < /esc < 
0 . 2 . 

Our mock constraints on ^Lw/^ion given 100 hours of in¬ 
tegration on a single sky region (EMI) can only rule out 
rather extreme cases. For example, this scenario rules out 
the 50,000 K toy stellar population with /esc ^ 0.2 at one 
extreme, and bottom-heavy IMFs with escape fractions be¬ 
low /esc ^0.1 at the other extreme. A stronger detection of 
turning point D, achieved by EM2, tightens these constraints 
considerably. The pure 50,000 K stellar population would 
require /esc ^ 0.01, while a stellar population with an over¬ 
abundance of lower mass stars would require /esc ^ 0.2. Note 
that the surface temperatures of PopIII stars are expected to 
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Figure 8. Left: Ratio of yields in the ionizing (/zv > 13.6 eV) and LW (10.2 < /?v/eV < 13.6) bands pe r stellar baryon as a f unction of metallicity and stellar 
IMF. Symbols represent model SEDs generated with STARBURST99 (those shown in Figures 1, 3, and 5 of ILeitherer et all 19991) . while the horizontal lines show 
the values one obtains for pure blackbodies at 10,000, 30,000, and 50,000 K from bottom to top. The filled symbols ai‘e investigated in more detail in the right 
panel. Right: Constraints on the stellar population and the escape fraction of ionizing radiation. The solid contour is the 2-a constraint on our reference model, 
i.e., identical to the green area of panel (e) in Figure[7] while dashed contours correspond to turning point constraints from EM2 (see Table^, which has a tighter 
constraint on the emission maximum (turning point D). The blue, green, and red bands have the same value of A^ion/M.w as the filled plot symbols in the left-hand 
panel, while the cross-hatched band instead adopts a pure 50,000 K blackbody spectrum for the stellar population. The width of each band corresponds to a 
factor of 2 change in the escape fraction, 0.1 < /esc ^ 0.2. 


be ^ 10^ K, which only strengthens our limits quoted for the 
50,000 K population. Our reference model assumes a typi¬ 
cal PopII stellar population, so it is reassuring to see that our 
constraints coincide with the blue diagonal band, which rep¬ 
resents a standard Salpeter IMF. 

Synthesis models for black hole populations are growing 
in maturity, though still only loosely cons trained by obser¬ 
vation s, especially at low metallicities (e.g.. lBelczvnski et al.l 
120081) . An immediate interpretation of will thus be very 
challenging barring progress on this front in the coming years, 
given that we cannot eliminate the degeneracy with the SFRD 
as we did previously by looking at the ratio ^Lw/^ion. rather 
than either quantity independently. For simplicity, we assume 
an a = 1.5 pow er-law spectrum abov e 0.2 keV consistent with 
the findings of iMineo et al.l ([20123), and /* =0.1. The 1-D 
marginalized PDFs for for EMI and EM2 are shown in 
EigurejOj Eactor of ^ a few enhancements are allo wed out 
to z < 4 (iDiikstra et al.ll201^ iBasu-Zvch et al.ll2013l) . though 
the redshifts probed by the global 21 -cm signal are far beyond 
the reach of the techniques used to establish such limits (the 
cosmic X-ray background and image stacking, respectively). 
All signal extraction models considered here rule out factor 
of 2 enhancements to fx at the ~ 3o level assuming /* =0.1. 
We will revisit this type of constraint in 0 

Our reference model is seemingly inconsistent with star for¬ 
mation in molecular halos and a stellar IME that yields more 
high-mass X-ray binaries than average (per unit star forma¬ 
tion). This does not rule out star formation in molecular halos 
or a top-heavy IME, it just rules out such sources as important 
drivers of the turning points. If we assume that PopIII stars 
have Alw = 4800, a SERD of « 3 x 10^^ Mq yr^' cMpc^^ 
would be required to match the constraint on provided by 
turning point B (following Eq. 17 of Paper I), which corre¬ 
sponds to ft fs 0.1 in Tjnin = 300 K halos. Such a population 
would have to die out rapidly in order for turning point C to 
be unaffected. Put another way, if PopIII stars do form rel¬ 
atively efficiently at z ~ 30, and continue to do so for more 



Figure 9. Constraints on sources of X-rays. Ttie solid and dastied liistograms 
are identical to tlie blaclc and green PDFs sliown in panel (i) of Figure [7] 
respectively, wliile dotted vertical lines denote values representative of 2x 
and 4x enhancements to the Lx-SFR relation, holding the star-formation ef¬ 
ficiency constant at /* =0.1. For reference, we show the 3 — 0 upper limit 
from extraction EM2 in red. Even if /* = 0.01, we limit fx < 20. 

than ^100 Myr, we should expect the position of turning 
point C to change (relative to our reference model) due to a 
stronger Ly-a background and potentially a stronger X-ray 
background, depending on the properties of PopIII remnants. 

4. DISCUSSION 

Our results suggest that simultaneous fits to the three spec¬ 
tral turning points of the global 21-cm signal can yield power¬ 
ful constraints on the properties of the Universe’s first galax¬ 
ies. A simple 4-parameter model can be constrained quite 
well in only 100 hours of integration on a single sky region, 
provided an ideal instrument. The t, parameters place inter¬ 
esting constraints on the properties of the first generations of 
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stars and black holes, while constraints on the characteristic 
redshift-dependent mass of star-forming galaxies follows im¬ 
mediately from constraints on Ty^m- In this section, we discuss 
these findings within a broader context, focusing in particu- 
la r on how our res ults depend on the assumed measurement 
( ^4. Il l and model ( ^4.21) . and how our fitting procedure might 
be improved to maximize the retu rn from ongoing and near¬ 
future observing campaigns ( ^4.3b . 

4.1. Are all three points necessary? 

Our forecasts have so far assumed that all three spectral fea¬ 
tures in the 40 < v/ MHz <120 window are detected and 
characterized reasonably well, apart from the EMI detection 
of turning point D which was only marginal. Given practical 
limitations in constructing an instrument with a broad-band 
response, the ionospheric challenges at low frequencies, and 
a weak emission feature potentially plagued by terrestrial ra¬ 
dio frequency interference (REI), it is worth asking: must we 
detect all three features at once to constrain even the simplest 
of galaxy formation models? 

Figures [To] and [TT] show the constraints on our 4-parameter 
model assuming only a subset of the turning points are de¬ 
tected. We consider all possible cases, except a scenario in 
which only turning point B is detected, as it seems unlikely 
that one could recover this feature from the foreground with¬ 
out help from neighboring spectral structure, given its ampli¬ 
tude of < 5 mK. Note that the black contours in each plot are 
identical to the 95% confidence regions in Figure |7] though 
the X and y scales of each individual panel here are much 
broader than those in Figure [T] due to the reduced quality of 
constraints. Blue contours denote fits including two turning 
points, while green cross-hatched regions correspond to fits 
including only a single turning point. Because the PDFs for 
the one- and two-point fits are broad, they tend to become 
noisy. This behavior is expected: by design, walkers spend 
less time in low-likelihood regions. If those regions of param¬ 
eter space are large (which they are for the one- and two-point 
fits), it will take a long time to properly explore them. 

In Figure [T0| we focus on the case in which the emission 
maximum, turning point D, is not used in the fit. In the most 
optimistic case, both turning points B and C are still detected, 
and give rise to the constraints shown in blue. As expected, 
constraints on are virtually nil except for a weak upper 
limit (panel g). However, constraints on ^lw, and Tjjiin 
remain largely intact. The subtle detours away from the black 
contours in panels b, f, and h toward small values of are 
real: they indicate scenarios in which heating is negligible 
and turning point C is induced by ionization (see §3.2.2 in 
Paper I). Such models would likely lead to an early end to 
the EoR and thus a large value of the Thomson optical depth, 
Tg, though without a detection of turning point D or a prior 
on Tg such scenarios remain allowed. In a more pessimistic 
scenario in which only the absorption minimum, turning point 
C, is detected, 2 — a constraints span ^ 3 orders of magnitude 
(green contours and cross-hatched regions), though still rule- 
out large regions of currently permitted parameter space. 

In the event that the lowest frequency feature, turning point 
B, is not detected, we instead arrive at the constraints shown 
in Figure [H] Provided that turning points C and D are both 
still detected, we obtain the blue contours, which are broader 
by ^ 1 order of magnitude in each dimension except for ^jon, 
though they still close within the broad space defined by our 
priors. If only the emission maximum is detected, we instead 
derive the contours in green. Tyni„ is unconstrained in this sce¬ 


nario, and only limits are available for each parameter when 
marginalizing over the others. 

The results shown in Figures [T0|and[TT]are almost certainly 
optimistic since it is the spectral structure which makes sig¬ 
nal extraction possible in the first place. With a narrow-band 
measurement of the global 21-cm signal that only includes 
two features, we should then expect the errors on the posi¬ 
tions of those turning points to be larger than those quoted 
in Table [T] Even so, such constraints would still be a big 
step forward, ruling out large regions of currently permitted 
parameter space and providing priors for next-generation ex¬ 
periments. 

4.2. Assumptions Underlying the Physical Model 

Our constraints on the four-parameter model of course as¬ 
sume that this model is “correct,” i.e., its parameters are as¬ 
sumed to be physically meaningful. In this section, we de¬ 
scribe the assumptions and approximations at the heart of this 
model and the circumstances in which they may deteriorate. 
This will provide a basis for our final piece of discussion re¬ 
garding the use of independe nt c onstraints and model selec¬ 
tion techniques, to follow in ^4.31 

4.2.1. The Star Formation History 

Our /coll-based recipe for the global 21-cm signal is cer¬ 
tainly not unique in its ability to model the first galaxies and 
the high-z IGM. For example, it would not be unreasonable 
to employ a mo re flexible “multi-population” model (e.g., 
lFurlanettol2006^ in which the Fy-a, soft UV, and X-ray back¬ 
grounds are produced by distinct sources, whose redshift evo¬ 
lution, photon production rates, and/or spectral energy distri¬ 
butions are allowed to be different. This approach may be 
warranted, given that the radiative properties and formation 
efficiencies (with time) of PopII and PopIII stars are expected 
to be different. 

Some authors have instead used empirical constraints on 
the SFRD at high-z to model the global 21-cm signal 
(lYaiima & KhochfMl2015h . While in principle such models 
are capable of more varied star-formation histories than our 
own, and can more seamlessly be compared to pre-existing 
empirical constraints on the SFRD in the post-EoR Universe 
(from which such SFRD models were first born), they have 
more free parameters and potentially obfuscate the dominant 
mode of star formation, which is of primary interest in this 
study. It would be straightforward to generalize our code to 
test empirically-calibrated parameterizations, which have the 
greatest strength at the lowest redshifts (z < 10), thus comple¬ 
menting the /coll approach, which is likely most accurate at 
the highest redshifts. 

Such changes to the underlying model would prevent some 
of the analysis so far presented. For example, our con¬ 
straints on the stellar IMF and escape fraction relied on the 
assumption of a single population well described by time- 
independent values of /*, the IMF (which we model implic¬ 
itly through A/w and Mon), and /esc. Such analyses could still 
be applied for a single-population model with an empirical 
SFRD, but for any kind of multi-population model. Equation 
[TSlno longer applies. In addition, the value of Mon/M,w may 
take on a new meaning, since it could probe A/w of PopIII 
stars that induce turning point B, and the Mon of more ordi¬ 
nary PopII star-forming galaxies responsible for driving turn¬ 
ing point D. 

4.2.2. Stellar Population Models 
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Figure 10. Constraints on the 4-parameter model, assuming the emission maximum is not detected. Solid black contours are the 95% confidence intervals of 
our reference fit to all three features. Blue contours are 95% confidence regions obtained when only turning points B and C ai'e used in the fit, while the green 
hatched regions show ares of parameter space excluded at 95% confidence if only turning point C used in the fit. When only upper or lower limits are available, 
we denote them with arrows in the marginalized 1-D PDFs. 


Even with perfect knowledge of the SFRD, properly in¬ 
terpreting Nion/N-Lw in terms of the stellar population re¬ 
quires robust predictions from synthesis codes, which aim 
to generate model SEDs as a function of time. Despite a 
long history and plenty of observational datasets to compare 
to, such codes are still being revised to account for updates 
in e.g., stellar atmosphere models ( e.g., line blanketing and 
NLTE effects Pauldrach et al.ll200lh and evolutionary tracks 
(lEkstrom et al. I2012t iGeorgy et al.ll2013l) . Such changes are 
pertinent to the reionization community, given their effects on 
the i onizing photon production efficiency per unit star forma¬ 
tion (jOitE^EiEHllloISISEElfli^^^SSlIlSill). However, 
further revision of stellar population synthesis models may 
not change how we model the global 21-cm signal, but rather 
change how we interpret constraints on the parameters of our 
mode l, particularly ^lw and ^ion. We will revisit this point in 

4.2.3. Stellar Remnants and X-ray Emission 

A complete stellar synthesis code would model the rem¬ 
nants of stars, in addition to stars themselves, if a compar¬ 
ison to datasets in the X-ray band were of interest. This is 
because neutron stars and black holes, when in binary sys¬ 
tems, are known to dominate the X-ray luminosity of star- 
forming galaxies (without active nuclei: iGrimm et £0120031: 


iGilfanov et al.l 120041: iMineo et al.l |20T23), while supernovae 
can provide yet another source of X-rays, either via inv erse 
Comp ton scattering off hot electrons in the remnant (lOhl 
IMl, or indirectly by mechanically heating the interstellar 
mediu m (ISM) which then emits thermal bremsstrahlung ra¬ 
diation (IMineo et al.ll2012bh . While we are not in the business 
of comparing model and measured X-ray spectra, we are con¬ 
cerned with modeling the X-ray SED of galaxies insofar as it 
affects the thermal history of the IGM. 

The modeling of compact object populations has become a 
growing industry in recent years, motivated in large part by 
the development of gravitational wave obse rvatories, contin¬ 
ued interest in ultra-luminous X-ray sources (iBelczvnski et al.l 
l2002t IBelczvnski et al.l^2008^ . and the likely ir nportance of 
comp a ct objects in reheating of the high-z IGM jPowe^et^ 
2009; iMirabel et al.l 1201 It Power et all 120131: iFragos et^ 
20131) . As in the case of stellar population modeling, the num¬ 
ber of compact objects and their mass distribution is expected 
to depend strongly on the metallicity. Unfortunately, the ob¬ 
servational data is sparse, especially at low metallicities, mak¬ 
ing it difficult to calibrate the models to local analogs of high-z 
galaxies. 

Whereas our forecast for the stellar IMF and escape frac¬ 
tion relied on the assumption of time-independent (but free to 
vary) values for /* and Tmin, our ability to constrain fx was 
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Figure 11. Constraints on the 4-parameter model, assuming turning point B is not detected. Solid black contours are the 95% confidence intervals of our 
reference fit to all three features. Blue contours are 95% confidence regions obtained when only turning points C and D are used in the fit, while the green hatched 
regions show ares of parameter space excluded at 95% confidence if only turning point D used in the fit. When only upper or lower limits are available, we denote 
them with arrows in the marginalized 1-D PDFs. 


intimately linked to the precise value of /*. Without more 
robust predictions for the X-ray yields of stellar populations, 
interpretation of fx will hinge on assumptions, or hopefully 
independent constraints, on the efficiency of star formation in 
high-z galaxies. Even if Lx-SER does not evolve much with 
redshift, uncertainties in the SED of X-r ay sources will cloud 
inferences drawn from the global 21-cm (lM^^ochal2014^ . Dis¬ 
entangling the normalization and spectral shape of the X- 
ray background will likely requi re independent measurements 
of the 21-cm powe r spectrum jPritchard & Eurlanettol 120071 : 
iPacucci et al.ll2014h . 

4.2.4. Cosmology and the Mass Function of Dark Matter Halos 

We have fixed cosmological parameters as well as parame¬ 
ters governing the halo mass functio n, adopting the most up- 
to-date values from Planck and the ISheth & TormenI (119991) 
form of the mass function throughout (computed using HMF- 
CALC: lMurrav et alJ2013h . Variations in the cosmological pa¬ 
rameters alone should be a secondary effect to all astrophys- 
ical processes we consider, but potentially discernible with 
observations of the dark ages (v ^ 20 MHz), prior to first- 
light. Variations in the cosmological parameters will also in¬ 
fluence the abundance of halos, though discrepancies in the 
halo mass function in the literature are known to be primarily 
due to differences in calibration of the fitting functions rather 


than u ncertainties in cosmological parameters (iMurrav et alJ 
120131) . at least at low redshifts (z < 2). Calibration of the mass 
function at high redshifts and for low-mass halos in which the 
first objects form is limited given the dynamic range needed 
to resolve small halos in large volumes. I f the mass func¬ 
tion at z > 20 deviates significantly from the lSheth & Torm^ 
(Il999h form, it would certainly affect the way we interpret 
Tjmn, and thus should be considered an important avenue for 
future work. 

4.2.5. The Two-Zone IGM Formalism 

Our entire procedure hinges on the ability to rapidly gener¬ 
ate model realizations of the global 21-cm signal, which has 
led us to a simple two-phase IGM formalism rather than more 
detailed (and expensive) numerical or semi-numerical simu¬ 
lations. Whereas simple models have been compared to nu- 
meric al simulations in the context of the 21-cm power spec¬ 
trum (IZahn et al.l201 Ih . and found to agree quite well, no such 
comparison has been conducted for global models. As far as 
we can tell, this is because there has yet to be a single nu¬ 
merical simulation capable of self-consistently generating a 
synthetic global 21-cm signal. Doing so will require high dy¬ 
namic range, capable of resolving the first star-forming halos, 
the radiation backgrounds they seed, in a domain large enough 
to be considered a global volume element. 
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Without a suite of numerical simulations to calibrate 
against, we have not attempted to attach any intrinsic uncer- 
tainty associated with ou r model, as was done recently by 
iGreig & MesingeJ (1201 5h in the context of the 21 -cm power 
spectrum. However, we do expect this formalism to be ac¬ 
curate over nearly the entire redshift range covered by our 
calculations (i.e., we do not use it solely out of computational 
necessity). The two-zone formalism operates best when HII 
regions are distinct and have sharp edges, and the heating and 
Ly-a is well-modeled by a uniform background. At turning 
point D, overlap between bubbles is likely minimal given that 
the volume filling factor of HII regions is small (Q ^ 0.2). In 
addition, their edges are likely sharp since fx is at most ^ a 
few, in which case X-ray binaries d o not enhance ioniza tion 
and heating by much on small scales (iKnevitt et al.ll2014l) . As 
a result, we do not have reason to suspect a breakdown in the 
formalism, at least for the reference model we have chosen. 

4.3. Priors and the Prospects for Model Selection 

Changes to the physical model, like those discussed in the 
previous section, generally fall into two categories; they either 
(1) change how we interpret the constraints on model param¬ 
eters of interest, or (2) fundamentally change the character¬ 
istics of the modeled signal. For example, improvements to 
synthesis models of stars and black holes will change how the 
^ parameters relate to the stellar IMF and properties of stellar 
remnants, and thus change how we interpret ^ values. But, 
so long as we still employ the four-parameter model, our con¬ 
straints on the values of ^ will not change. If instead we intro¬ 
duced new parameters that allowed ^ or Tmin to evolve with 
redshift, we have then enhanced the flexibility of the model 
enough that we may now be capable of generating realiza¬ 
tions of the global 21-cm that our previous approach simply 
could not have. 

A “double reionization” scenario, which could lead to two 
emission featur es rather than our single “ turning point D,” is 
an implausible (iFurlanetto & LoebI120051) but illustrative ex¬ 
ample to consider in this context. Our four-parameter model 
simply could not produce two emission features. One could 
imagine less drastic changes that might still have new and po¬ 
tentially discernible effects on the signal through modulations 
of its shape, such as redshift-dependent ^ and r^in, feedback, 
and/or multiple distinct source populations. 

We should expect that more complex model parameteriza- 
tions will only have an easier time fitting the turning points, 
and thus a fit to the turning points alone may not enable one 
to constrain additional parameters. Use of a more complex 
parameterization may still be warranted if independent con¬ 
straints on one or more of the model parameters are available 
to be used as priors in the fit. However, if we do a “single- 
stage” fit, in which we fit a physical model directly to the data 
rather than using a computationally inexpensive intermediary 
to extract the turning points, we may find that a more complex 
model is required by the data. In order to justify the additional 
parameters rigorously, more advanced i nference tools are 
required (e.g., mul tinest,PolyChord; iFeroz et all 120091; 
iHandlev et al.ll20l3h to compute the Bayesian evidence. 

To date, there has only been one pape r on model selection 
for global 21-cm datasets (lHarkeril2 015) . The evidence is ex¬ 
pensive enough to compute that lHarkCT (120151) was limited to 
relatively low-dimensional spaces and simplistic signal mod¬ 
els. In the future, such tests will be required in order to test 
whether or not more complex models (i.e., those more com¬ 
plex than our four-parameter reference model) are required 


by the data. This presents a unique and challenging problem 
for ongoing and upcoming experiments and their associated 
signal extraction pipelines. 

5. CONCLUSIONS 

This work represents the first attempt to forecast constraints 
on astrophysical parameters of interest from mock observa¬ 
tions of the global 21-cm signal. There is clearly much still to 
be learned, even from synthetic datasets, about how observa¬ 
tions in (v, 87),) space translate to constraints on the properties 
of the IGM and/or the properties of high-z galaxies. Assum¬ 
ing an idealized instrument, signal recovery consistent with 
the values quoted in Table [T] and the validity of our four- 
parameter model for the global 21-cm signal, we find that; 

1. Constraints on the turning points constrain the model 
parameters well (to ^0.1 dex each, on average), with 
factor of ~ 2 improvements within reach of exper¬ 
iments viewing multiple sky regions and employing 
more complex foreground parameterizations (Figure|7]i. 
Such measurements would simultaneously constrain 
the ionization and thermal state of the IGM (Figures 
I5ll6t . perhaps most notably providing strong evidence 
for the beginning of the EoR at z 12. 

2. Our fiducial realization of the signal is inconsistent with 
star-formation in halos with virial temperatures below 
~ 10^’^ K at the 2-a level for the most pessimistic signal 
extraction scenario we consider. Such constraints are 
enabled in large part by a broad-band measurement of 
the signal, since affects all three turning points in 
the ^ 40 — 120 MHz interval (Figures [T]|2] and|7]l. 

3. In the simplest case, in which all model parameters are 
assumed to be constant in time, we can provide limits 
on both the escape fraction and the stellar IMF, primar¬ 
ily ruling out scenarios in which UV photons originate 
in extreme environments with very top-heavy IMFs or 
very high escape fractions (Figure|8ll. 

4. Our constraints on X-ray sources are comparable to 
those achieved at z < 4 via stacking and the cosmic X- 
ray background, though at z 20 to which the afore¬ 
mentioned techniques are insensitive (Figure |9]l. In the 
absence of independent information, this constraint re¬ 
quires an assumption about the star formation efficiency 
and X-ray SED of galaxies, however. 

5. With only a subset of the turning points, constraints 
on our reference model are considerably worse (Eig- 
uresfTOlandfTTTi. The lowest frequency features (turning 
points B and C) hold the most power to constrain Zinin, 
which will make it difficult to constrain and ^ion 
with observations confined to the highest frequencies. 
Isolated detection of the absorption feature is the most 
valuable single-point measurement, as it leads to confi¬ 
dence contours which close over the prior space, except 
in the case of ^ion. 
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